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With the help of von Neumann entropy, we study numerically the localization properties for two 
interacting particles (TIP) with the on-site interaction in one-dimensional disordered, quasiperiodic 
and slowly varying potential systems, respectively. We find that for TIP in disordered and slowly 
varying potential systems, the spectrum-averaged von Neunmann entropy (E v ) first increases with 
interaction U until its peaks, then decreases as U gets larger. For TIP in the Harper model, 
the functions of {E v ) versus U are different for particles in extended and localized regimes. Our 
numerical results indicate that for these two-particle systems, the von Neumann entropy is a suitable 
quantity to characterize the localization properties of particle states. Moreover, our studies propose 
a consistent interpretation of the discrepancies between previous numerical results. 
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I. INTRODUCTION 

Interacting uniform electronic systems^ and non- 
interacting disordered electronic systems 2 ^ are two of the 
most intriguing, albeit difficult, subjects in condensed 
matter physics. As a result of the complexity of the 
simultaneous presence of randomness and interactions, 
few definitive results are known4£ To understand the ef- 
fects of electrons interaction on the localization proper- 
ties in a random potential, Shepelyansky proposed some 
years earlier that it would be worthwhile to consider 
the simple case of two interacting particles(TIP) in a 
one-dimensional (ID) random potential. 6 Subsequently, 
extensive efforts have been devoted to TIP in various 
systems ^ i 8 ' 9 : 10 ! 11 ! 12 ! 14 In particular, two interacting elec- 
trons with total spin zero in a ID Harper model,— in 
Fibonacci and Thue-Morse lattices^ are studied. The 
behavior of TIP has been studied using the time evo- 
lution of wave packets ^ii 8 .^ exact diagonalization j2iifi 
Green functio n 11 : 12 ! 13 and transfer-matrixjii etc. Due 
to different definitions of the localization length for TIP 
and different methods applied, there are discrepancies 
about the effect of the combination of disorder and in- 
teraction. For example, for TIP in ID disordered po- 
tentials, Shepelyansky- 6 - and other a 11 ' 12 have found that 
for small interaction strengths, the two-particle interac- 
tion can enhance the Anderson localization length, while 
Romer and Schreiber 14 found no enhancement of the lo- 
calization length when the system size grows to infin- 
ity. Evangelou et al. have pointed out that stronger lo- 
calization occurs at large interaction when compared to 
the noninteracting caseJ^ At the same time, for TIP in 
a ID Harper model, it was found that the interaction 
would induce localization effect? 7 - : 15 : 16 while Evangelou 
and Katsanos^ found that the effects of particle interac- 
tions are different for electrons in extended and localized 
regimes. 



On the other hand, quantum entanglement, which 
attracting much attention in quantum information, 17 
has been extensively applied in condensed matter 
physicsj 18 ' 19 : 20 ' 21 : 22 : 23 For examples, quantum entangle- 
ment measured by the von Neumann entropy was stud- 
ied in the Hubbard model for the dimer cas&i^, in the 
extended Hubbard model for different band filling^, in 
quantum small- world networks^ 2 -, and in low-dimensional 
semiconductor systems^ 3 -. It was found that the von 
Neumann entropy is suitable for analyzing the inter- 
play between itinerant and localized features^ 8 -, as well 
as characterizing quantum phase transitio n 20 : 21 and the 
localization-delocalization transition of electron states^. 

In this paper, we perform a detailed study of the 
von Neumann entropy for TIP in ID disordered and 
quasiperiodic systems respectively, taking into account 
different on-site interactions U at various on-site poten- 
tial strengths. We find that for TIP in disordered and 
slowly varying potential systems, the spectrum-averaged 
von Neunmann entropy (E v ) first increases with interac- 
tion U until its peaks, then decreases as U gets larger. 
For TIP in the ID Harper model, the functions of (E v ) 
versus U are different for particles in extended and lo- 
calized regimes. Finally we study a two-particle system 
based on the slowly varying potential model i 24 : 25 : 26 From 
these studies, we can conclude that for TIP systems, the 
von Neumann entropy is a suitable quantity to charac- 
terize the localization properties of particle states. 

The paper is organized as follows. In the next section 
the von Neumann entropy is described. In Sec. IIIII the 
numerical results for TIP in ID disordered, quasiperi- 
odic, and slowly varying potential systems are presented, 
respectively. And we present our conclusions and discus- 
sions in Section IIVI 
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II. TIP MODELS AND VON NEUMANN 
ENTROPY 

A. TIP models 

Following recent literature 6 -^ ' 12 ' 14 ' 15 ' 16 the eigenvalue 
equation for the TIP in ID system can be written as 

(^ni ~t~ ^^TOl ,712 /V^ll ,ri2 ^(^711 + 1,712 

)=El/j nu n 2 , (1) 

where e n is the on-site potential, t is a nearest-neighbor 
hopping integral and £/ characterizes the on-site interac- 
tion between particles. 

Eq. U) actually can describe the behaviors both bosons 
and fermions, i.e., two spinless bosons or two electrons 
with opposite spins. 

For spinless bosons with the on-site interaction, the 
tight-binding Hamiltonian can generally be described by 

JV-l N 
H = * Y] (4iCn+l + C f n+1 C n ) + V E n C n C n 



n=l 



n=l 



N 



(2) 



where & n (c n ) is the boson creation(annihilation) opera- 
tor of the nth site. The generic eigenstate for two spinless 
bosons is the superposition 



N 



E c 



ni<ri 2 



N 

„ 2 K,n 2 ) = E d,T. a 4i c L 1°) ' 

ni<n 2 



(3) 

where |0) is the vacuum and „ 2 is the amplitude of 
wave function. From Eqs.Q and |5]) we can obtain the 
eigenvalue equation (fTJ). 

For electrons with the on-site interaction, the tight- 
binding Hamiltonian can be described b y 8 ' 9 ' 10 

N 

H = t E ( C li,cr c n+l,a- + <4i+l,CT C n,<0 
n—1 (7 

N N 
+ EE 

(4) 



where c^ n a {c n ,a) is the electron creation(annihilation) 
operator for the electron at the nth site with spin a = 
±i. For two electrons the Hilbert space can be conve- 
niently divided into a singlet subspace with total spin 
5 = and a triplet subspace with total spin S = 1, re- 
spectively. Since the triplet subspace permits no double 
occupation, it is not affected by the on-site interaction. 
In order to analyze the effect of the on-site interaction 
U, we will consider the case of the two electrons with 
opposite spins, i.e., the singlet subspace. In a chain with 
N sites the singlet subspace is spanned by N(N + l)/2 



spatially symmetric basis functions 

_ / 75( c Lt c Ll + c l 2 ,Ax,0 1°) for ni < ™ 2 ' 
- „t J 



( c ni,T C n 3 ,,l) 1°) for "1 = "2, 



(5) 



where |0) is the vacuum. An eigenstate for two electrons 
with the spatially symmetric wave functions is in general 
the superposition 



A' 



= E ^.na 



n 2 ) 



(6) 



ni <n 2 



where <f>" n is the amplitude of wave function. By mak- 
ing the transformation 



ib a 



V^^ni,n a for n l = n 2, 



the eigenvalue equation obtained from Eqs. (|H— [5]) can 
be written as equation (JTJ) . 

In the following our numerical method is described for 
bosons. The extension for fermions is straightforward. 



B. von Neumann entropy 

For the two particles in the system we are studying, 
there are three local states at each site, \2) n , |l) n , |0) n , 
corresponding to the state with two, one or zero particles 
at the nth site, respectively. The local density matrix p n 
is define d 18 ' 20 by 

Pn = Z2n I 2 ) nn (2| + Z\ n |1) „„ (1| 

+ (1 - z ln - z 2n ) |0) „„ <0| . (7) 
For two spinless bosons, 



Z 2n = (a\ clcnclcn \a) = 



n . n 



(8) 



and 



zi n = (a\ c n c n | a) - 2z 2n 

N N 

= E <m<*m + E C',„«,n- (9) 
m(m>n) m' (m' <n) 

Consequently, the corresponding von Neumann entropy 
related to nth site is 

E vn = "(I - z ln - Z 2n ) fog 2 (l - Zln ~ Z 2n ) 

-Zi n log 2 Zi„ - z 2n \og 2 z 2n . (10) 

For nonuniform systems, the value of E" n depends on 
the site position n. At an eigenstate a, we define a site- 
averaged von Neumann entropy 



1 N 

n=l 



(11) 
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For a delocalized state that all ib" _ = 1 for 

rni '" 2 y/N(N+l)/2 

all m < ri2, this definition gives E£ « ^ log 2 at large 
JV, while for a localized state that ^>"o„o = 1 at given 
and n%, E™ = 0. In this paper, all the values of E" are 
scaled by ■£> log 2 From the two examples, we know 
that the scaled E" is close to 1 for eigenstates are ex- 
tended and almost vanishes for eigenstates are localized. 
Henceforth, we omit "scaled" for simplicity. 

In order to analyze the influence of system parame- 
ters like the on-site interaction U, on the von Neumann 
entropy for all the eigenstates, we define a spectrum- 
averaged von Neumann entropy as a further gross mea- 
sure 

OL 

where M is the number of all the eigenstates. 



III. NUMERICAL RESULTS 

From now we consider only the repulsive interaction 
(U > 0). We directly diagonalize the eigenvalue Eq.|T|) 
with the periodic boundary condition at finite system 
sizes and obtain all eigenvalues E a and the corresponding 
eigenvectors |a). Without loss of generality, the hopping 
integral t is taken as units of energy. From the formulas 
(f7l — TT2|) . we then can obtain the site-averaged von Neu- 
mann entropy E" and the spectrum-averaged von Neu- 
mann entropy (E v ), respectively. 



A. TIP in a disordered potential chain 

For TIP in a disordered potential chain, the on-site 
potential e n in Eq.Q are random variables uniformly 
distributed with \—W, W]. Here W characterizes the de- 
gree of on-site disorder as in the Anderson models For 
this model, we calculate the spectrum-averaged von Neu- 
mann entropy (E v ) with a given set of parameters W and 
U. For every set of parameters W and U, the disorder 
average is taken over 100 samples*^ More samples simply 
give similar results. 

For ID Anderson model in the absence of the 
interaction([7 = 0), it is well known that all the eigen- 
states are localized and the one-particle localization 
length is £ w 25t 2 /W 2 at the energy band center^ 
FigH] gives the spectrum-averaged von Neumann entropy 
(E v ) versus the disorder parameter W . It shows that 
(E v ) monotonically decreases as W increases, reflect- 
ing the trivial localization effect of the on-site disor- 
der in the model. To find the correlation between (E v ) 
and the localization properties for TIP systems, we also 
study the inverse participation ratio (IPR)^, defined by 

JV 

= ( N{N 2 +1) J2 Id.nJ 4 )^ wnicn S ives the ratio of 

m<ri2 
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FIG. 1: The spectrum-averaged von Neumann entropy 
(E v ) versus the disorder W for U = at JV = 90. 




FIG. 2: The relation between the site-averaged von 
Neumann entropy E" and the corresponding £ Q at 
W = 1 and U — for a typical sample. The E£ versus 
log 10 £ a is plotted in the inset. 



lattice sites occupied by particles to all lattice sites at 
an eigenstate a. The larger £ Q is, the more delocalized 
the eigenstate is. In Fig[2] we plot the site-averaged von 
Neumann entropy E™ versus £ Q at U = for a typical 
sample with W = 1. On the whole, E" increases loga- 
rithmically with as can been seen in the inset of Figf21 
so von Neumann entropy can well reflect the localization 
properties of two-particle eigenstates. 

For U > 0, the spectrum-averaged von Ncunmann en- 
tropy (E v ) as a function of U at W — 1, 2, 3 are plotted 
in FigO The results are similar for other W . For all W, 
we find that (E v ) first increases until its peak as U in- 
creases from zero, then decreases as U gets larger. There 
is a U*, at which (E v ) is equal to that at U = 0. When U 
is smaller(greater) than U* , (E v ) is larger (smaller) than 
the value of (E v ) at U = 0. To understand the effect of 
the interaction U on the von Neumann entropy, we cal- 
culate the site-averaged von Neumann entropy E" and 
illustrate the results in Figf4]at different U for W = 1. It 
shows that only a small portion of the two-particle eigen- 
states have its value of E™ changed by the interaction U. 
For small U, most of the newly created eigenstates are 
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FIG. 3: The spectrum-averaged von Neumann entropy 

(E v ) as a function of interaction strengths U for 
(a,)W = l,(b)W = 2 and (c)W = 3, respectively. Here 
N = 90. 



in the main band and the corresponding E" are larger 
than that in the noninteracting case, while for large U, 
most of them are above the top of the main band and 
the corresponding E™ are smaller than that for particles 
without interaction. Therefore, the varying (E v ) with U 
is not monotonic as shown in Fig[3] These two competing 
effects reach an equilibrium at U = U* . 

For TIP in ID disordered potentials with U < t 
Shepelyansky 6 - and other o 11 ' 12 have found that the An- 
derson localization length increases with U . It agrees 
with our conclusions that for small U, (E v ) increases 
with U. For large interaction U, Evangelou et al^- have 
pointed out that stronger localization occurs when com- 
pared to the noninteracting case, because at U » t in- 
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FIG. 4: The site-averaged von Neumann entropy E™ 
and the corresponding eigenenergy E a at different U for 
a typical sample at W — 1. Here N — 90 and the 
spectrum-averaged von Neumann entropy 
(E v ) = 0.8431, 0.8556 and 0.7916 for U = 0, 2 and 6, 
respectively. 



teraction can significantly modify the energy spectrum 6 
and the created eigenstates have the comparatively small 
localization lengths^ The result is consistent with ours 
that for U > U* , the values of (E v ) are smaller than that 
at U = 0. 



B. TIP in a quasiperiodic potential chain 

From Eq.([T]), the eigenvalue equation for TIP in a 
quasiperiodic potential chain based on the Harper model 
can be described by££ 

[Acos(27rcrni -I- f3\) + A cos(27rern 2 + (3 2 ) + U8 ni , n2 ]ip ni}n2 

+t(lp ni +l,n 2 + 1pni-l,n 2 + ^ni,n 2 +l + VVi.na— l) 
= Etp nuri2) (13) 

here the parameter A characterizes the strength of the 
quasiperiodic potential, a and /3\,2 are constants. As a 
typical case, we set a = — l)/2, /3i = fa = 0. As 
is customary in the context of quasiperiodic system, the 
value of a may in fact be approximated by the ratio of 
successive Fibonacci numbers — F n — F„_2 + F n -i . In 
this way, choosing a — F n -\/F n ~ (\/5 — l)/2 and the 
system size N = F ni we can obtain the periodic approxi- 
mant for the quasiperiodic potential. In our calculation, 
N is chosen as Fibonacci numbers 34, 55 and 89, respec- 
tively. 

In the absence of the interaction U it is foun d 28 ' 29 that 
that for A < 2 the spectrum becomes continues and all 
eigenstates are extended. For A > 2 the spectrum is 
pure point and all eigenstates are exponentially local- 
ized. For A = 2 the situation is critical with a singular- 
continuous multifractal spectrum and power law local- 
ized eigenstates. Metal-insulator transition can occur at 
A = 2. Obviously they are different from these for the 



5 



0.8 



A 

iu > 0.6 

V 



0.4 



0.2, 



0.98 



11=0 



o 

o 
o 
o 
o 



"0 0.5 1 1.5 



2.5 3 3.5 4 



(a) 




(b) 

FIG. 5: (a)The spectrum-averaged von Neumann 
entropy (E v ) as a function of quasiperiodic potential 
strength A at U = 0. Here N = 89. (b)The (E v ) versus 
A at U = for different N. 



ID Anderson model where all eigenstates are localized 
for W ^ 0. 

Figf5)Ja) shows the spectrum-averaged von Neumann 
entropy (E v ) with respect to the quasiperiodic potential 
strength A at U = 0. We observe that there is a sharp 
decrease of (E v ) near the critical value A c = 2. For A < 
A c which corresponds to extended states, all (E v ) are 
around 1, while for A > A c corresponding to localized 
states, all (E v ) are far less than 1. In FigfSJb) we plot 
(E v ) as a function of A at U = for N = 34, 55 and 89. 
The curves cross at A « A c . The crossing point separates 
the extended(A < A c ) and the localized(A > A c )regimes. 

For U > 0, the spectrum-averaged von Neumann en- 
tropy (E v ) as a function of U at different A are shown 
in Fig|5J Here we choose A = 1, 2 and 3 as examples for 
extended, critical and localized regimes, respectively. At 
A = 1 and A = 2, (E v ) monotonically decreases when U 
increases from zero. At A = 3 (E v ) increases for small U 
and decreases for large U, respectively, which is similar 
to that shown in Figj3] for TIP in a disordered potential 
chain. For A in extended and localized regimes, (E v ) as 
a function of U are similar to that for A = 1 and A = 3, 
respectively. We also find that at the extended and crit- 
ical regimes, the newly created eigenstates due to the 
interaction have small E£ , which will reduce the value of 
(E v ). At the localized regime, the effect of interaction 
U on (E v ) is similar to that in the disordered potential 
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FIG. 6: The spectrum-averaged von Neumann entropy 
(E v ) as a function of interaction strengths U for 
(a)A = l,(b)A = 2 and (c)A = 3, respectively. Here 
N = 89. 



system, so (E v ) increases for small U, decreases for large 
U, respectively. 

For TIP in the ID Harper model, Shepelyansky et 
a? ] 7 i 15 i 16 found that the interaction would induce local- 
ization effect for all A . At the same time, Evangelou 
and Katsanos^ found that, in the extended(A < 2) and 
the critical(A = 2) regimes, the velocity and the diffu- 
sion coefficient of TIP will decreases due to the localized 
pairing states. In the localization regime(A > 2), they 
found that propagation enhancement for small interac- 
tion and strong localization for large interaction, as in 
disorder systems. Apparently there are discrepancies in 
their results at the localization regime. Comparing to 
our results, in the extended and the critical regimes, we 
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FIG. 7: The spectrum-averaged von Neumann entropy 

(E v ) as a function of potential strength A at U = 0. 
The dashed line is linearly fitted for the corresponding 
data at < A < 2. Here ira = 0.2, v = 0.7 and N = 89. 

find (E v ) decreases with U, which agrees with both of 
their conclusions that the interaction can induce localiza- 
tion effect 15 ' 16 In the localized regime, we find (E v ) 
increases for small U and decreases for large U, respec- 
tively which is consistent with the results of Evangelou 
and Katsanos^ that the interaction has different effects 
on localization properties at small and large U. 

C. TIP in a slowly varying potential chain 

Nest we study TIP moving in a ID system based on 
the slowly varying potential model , 24 i 25 i 26 From Eq.fTT]), 
the eigenvalue equation can be described by 

[Acos^cm^ + Pi) + A cos(7ran2 + P2) + U5 ni . n2 ]ipm.n 2 

+t(fp ni +l,n 2 + - 0ni-l,n 2 + VVii,n 2 + l + V , ni,n 2 -l) 

= Et/i nun2 , (14) 

here A, a and v are positive numbers. For a irrational 
with v > 2 or v = 1, this is equivalent to the models dis- 
cussed in Sec. Mil (A) and (B), respectively. For < v < 
1, in the absence of the interaction U, it is well known 24 
that there are two mobility edges at E c = ±(2.0 — A) 
provided A < 2.0. It is found that extended states are in 
the middle of the band ( \E\ < 2.0 — A ) and localized 
states are at the band edge ( 2.0 - A < \E\ < 2.0 + A ). 
In other words, for A < 2.0, the extended and localized 
eigenstates coexist in contrast to the models studied in 
Sec.|TTIl(A) and (B). For A > 2.0, all states are found to 
be localized. So in this model there are always localized 
eigenstates at A ^ 0. Obviously the spectrum properties 
are different from that for the Anderson model and the 
Harper model. 

Fig[7| shows the spectrum-averaged von Neumann en- 
tropy (E v ) as a function of A at U = 0. For < A < 2 
the data can be well fitted into a line. For A > 2 all 
data points lie far away from the proposed line, i.e., there 
is an abrupt decrease in (E v ) at the critical parameter 



FIG. 8: The spectrum-averaged von Neumann entropy 
(E v ) as functions of interaction strengths U at different 
A. Here na = 0.2, v = 0.7 and N = 89. 



A = 2. The critical parameter separates the localized 
regime(A > 2) and the regime where the localized and 
extended eigenstates coexist(A < 2). This result is con- 
sistent with the spectrum properties for the model, 24 so 
the spectrum- averaged von Neumann entropy (E v ) can 
describe the localization properties for two-particle sys- 
tems. 

Fig[8] shows the spectrum-averaged von Neumann en- 
tropy (E v ) versus the interaction U at A < 2. The re- 
sults are similar to that for A > 2. We again find that 
for all A, as U increases, (E v ) increases for small U, de- 
creases for large U, respectively, similar to that shown in 
Fig|3]for TIP in a disordered potential chain. It indicates 
that propagation enhancement for small interaction and 
strong localization for large interaction. In other words, 
as long as there are localized eigenstate at single-particle 
models, the effect of interaction on the localization prop- 
erties of TIP is similar to that in disorder systems. 



IV. CONCLUSIONS AND DISCUSSIONS 

With the help of von Neumann entropy, we have stud- 
ied the effect of the on-site interaction U on the local- 
ization properties of TIP in ID disordered, quasiperiodic 
and slowly varying potential systems, respectively. 

For TIP in a disordered potential chain and a slowly 
varying potential chain , as U increases from zero, we find 
that at first (E v ) increases for small U, then decreases 
as U gets large. It means that there are propagation 
enhancement for small interaction and stronge localiza- 
tion for large interaction. For TIP in ID Harper model 
with particles in the extended(A < 2) and critical(A = 
2) regimes, (E v ) decreases as U increases, which indi- 
cates that the interaction would induce localization ef- 
fect, while in the localized regime(A > 2), the (E v ) as a 
function of U is similar to that for TIP in the disordered 
potential chain. From our studies, we find that (E v ) is a 
suitable quantity to describe the localization properties 
for two-particle systems. 
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Summarizing all results from the three models, we can 
conclude that provided localized cigenstates at single- 
particle case exist, the delocaliation (localization) effect 
happens for small interactions (large interactions), while 
single particle states are extended or critical, the inter- 
action always induces the localization effect. According 
to our results, we propose a consistent interpretation of 
the discrepancies between previous numerical results. 



ful discussions. This work is partly supported by Na- 
tional Nature Science Foundation of China under Grants 
Nos. 90203009, 10175035, and 10674072, by SRFDP 
No.20060319007, by EYTP of MOE, China, by NSF of 
Jiangsu No. 06KJD140135, and by the Foundation of 
NJUPT No. NY205050. 



Acknowledgments 

We would like to thank Zhongxia Yang and Shijie 
Xiong for a critical reading of the manuscript and use- 



* Email address:lygong@njupt. edu.cn. 

t Corresponding author. Email ad- 

dress:pqtong@njnu. edu.cn. 

1 E. H. Lieb and F. Y. Wu, Phys. Rev. Lett. 20, 1445 (1968). 

2 P. W. Anderson, Phys. Rev.109, 1492 (1958). 

3 B. Kramer and A. MacKinnon, Rep. Prog. Phys. 56, 1469 

(1993) . 

4 P. A. Lee and T. V. Ramakrishnan, Rev. Mod. Phys. 57, 
287 (1985). 

5 D. Belitz and T. R. Kirkpatrick, Rev. Mod. Phys. 66, 261 

(1994) . 

6 D. L. Shepelyansky, Phys. Rev. Lett. 73, 2607 (1994). 

7 D. L. Shepelyansky, Phys. Rev. B 54, 14896 (1996). 

8 S. N. Evangelou, D. E. Katsanos, Phys. Rev. B 56, 12797 
(1997). 

9 P. E. de Brito, E. S. Rodrigues,and H. N. Nazareno, Phys. 
Rev. B 73, 014301 (2006). 

10 S. N. Evangelou, Shi-Jie Xiong, and E.N. Economou, Phys. 
Rev. B 54, 8469 (1996). 

11 F. von Oppen, T. Wettig, and J. Muller, Phys. Rev. Lett. 
76, 491 (1996). 

12 P. H. Song and Doochul Kim, Phys. Rev. B 56, 12217 
(1997). 

13 K. Frahm, A. Miiller-Groeling, J.-L. Pichard, and D. Wein- 
mann, Phys. Rev. Lett. 78, 4889(1997). 

14 R. A. Romer and M. Schreiber, Phys. Rev. Lett. 78, 515 
(1997). 

15 A. Barelli, J. Bellissard, P. Jacquod and D. L. Shepelyan- 
sky, Phys. Rev. Lett. 77, 4752(1996). 



16 A. Barelli, J. Bellissard, P. Jacquod and D. L. Shepelyan- 
sky, Phys. Rev. B 55, 9524 (1997). 

17 See, for example, The Physics of Quantum Information, 
edited by D. Bouwmeester, A. Ekert, and A. Zeilinger 
(Springer, Berlin, 2000). 

18 P. Zanardi, Phys. Rev. A 65, 042101 (2002). 

19 Yu Shi, J. Phys. A 37, 6807 (2004). 

20 S.-J. Gu, S.-S. Deng, Y.-Q. Li, and H.-Q. Lin, Phys. Rev. 
Lett. 93, 086402 (2004). 

21 D. Larsson and H. Johannesson, Phys. Rev. A 73, 042320 
(2006). 

22 L. Y. Gong and P. Q. Tong, Phys. Rev. E 74, 056103 
(2006). 

23 F. Buscemi, P. Bordone, and A. Bertoni, phys. Rev. A 73, 
052312 (2006). 

24 S. Das Sarma, S. He, and X. C. Xie, Phys. Rev. Lett. 61, 
2144 (1988). 

25 M. Griniasty and S. Fishman, Phys. Rev. Lett. 60, 1334 
(1988). 

26 D. J. Thouless, Phys. Rev. Lett. 61, 2141 (1988). 

27 For random potentials, we use the following definition of 
the spectrum-averaged von Neumann entropy: (E v ) = 
k 

JiJ2 E v = -r J2( jiJ2 E v )k, here X denotes the disorder 

a k — 1 a 

average, and K the number samples. 

28 S. Aubry and G. Andre, Ann. Isr. Phys. Soc. 3, 133 (1988). 

29 J. B. Sokoloff, Phys. Rep. 126, 189 (1985). 



APS/123-QED 



Von Neumann entropy and localization properties of two 
interacting particles in one-dimensional nonuniform systems 

Longyan Gong 1 > 2 > C * and Peiqing Tong 2 Q 
1 Department of Mathematics, Nanjing Normal University, 
Nanjing, Jiangsu, 210097, P. R.China 
2 Department of Physics, Nanjing Normal University, 
Nanjing, Jiangsu, 210097, P. R.China 
3 Department of Mathematics and Physics, 
Nanjing University of Posts and Telecommunications, 
Nanjing, Jiangsu 210003, P.R.China 
(Dated: today) 

Abstract 
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I. INTRODUCTION 



Inter acting uniform electronic systems! if] and non-interacting disordered electronic 

n n 

systems |2J, |3|] are two of the most intriguing, albeit difficult, subjects in condensed mat- 
ter physics. As a result of the complexity of the simultaneous presence of randomness and 
interactions, few definitive results are known. To understand the effects of electrons 
interaction on the localization properties in a random potential, Shepelyansky proposed 
some years earlier that it would be worthwhile to consider the simple case of two interact- 



ing particles(TIP) in a one-dimensional (ID) random potential. 



efforts have been devoted to TIP in various systems. 



10 



6 1 Subsequently, extensive 



11 



two interacting electrons with total spin zero in a ID Harper model, |8J] in Fibonacci and 



f 



14j In particular, 



Thue-Morse lattices 191 are studied. The behavior of TIP 



evolution of wave packets 



7 



ras been studied using the time 



8|, |l0j Green function 



11 



,112 



, |9j exact diagonalization 
and transfer- matrix. |l4j] etc. Due to different definitions of the localization length for TIP 
and different methods applied, there are discrepancies about the effect of the combination of 
disorder and interaction. For example, for TIP in ID disordered potentials, Shepelyansky 
and others 



11 



12| have found that for small interaction strengths, the two-particle interac- 



tion can enhance the Anderson localization length, while Romer and Schreiber 14| found no 
enhancement of the localization length when the system size grows to infinity. Evangelou et 
al. have pointed out that stronger localization occurs at large interaction when compared to 
the noninteracting case. 1 101 At the same time, for TIP in a ID Harper model, it was found 

ann 

that the interaction would induce localization effect, p], [15|, [16( while Evangelou and Kat- 
sanos [8| found that the effects of particle interactions are different for electrons in extended 
and localized regimes. 

On the other hand, quantum entanglement, which attracting much attention in quantum 



in 



22 



ormation. jl7] has been extensively applied in condensed matter physics. 



18 



19 



20. 



21 



231 ] For examples, quantum entanglement measured by the von Neumann entropy was 



studied in the Hubbard model for the dimer case [18j, in the extended Hubbard model 
for different band filling 2(|, in quantum small- world networks 22]. and in low-dimensional 
semiconductor systems 23J| . It was found that the von Neumann entropy is suitable for 



analyzing the interplay between itinerant and localized features 



181 ] , as well as characterizing 



quantum phase transition [20j, l21[ and the localization-delocalization transition of electron 



2 



states [22]. 

In this paper, we perform a detailed study of the von Neumann entropy for TIP in 
ID disordered and quasiperiodic systems respectively, taking into account different on-site 
interactions U at various on-site potential strengths. We find that for TIP in disordered 
and slowly varying potential systems, the spectrum-averaged von Neunmann entropy (E v ) 
first increases with interaction U until its peaks, then decreases as U gets larger. For TIP in 
the ID Harper model, the functions of (E v ) versus U are different for particles in extended 



and localized regimes. 



potential model . 



24 



25 



inally we study a two-particle system based on the slowly varying 



261 ] From these studies, we can conclude that for TIP systems, the 



von Neumann entropy is a suitable quantity to characterize the localization properties of 
particle states. 

The paper is organized as follows. In the next section the von Neumann entropy is 
described. In Sec. [TTT] the numerical results for TIP in ID disordered, quasiperiodic, and 
slowly varying potential systems are presented, respectively. And we present our conclusions 
and discussions in Section [TV] 

II. TIP MODELS AND VON NEUMANN ENTROPY 



A. TIP models 



Following recent literature 0, Q, Ijl, 
system can be written as 



14 



15 



161 ] the eigenvalue equation for the TIP in ID 



(—711 ^n>2 USni,n2)lpni,ri2 ^(V , ni+l,ri2 
-H0nx-l,Ti2 



(1) 



where e n is the on-site potential, t is a nearest-neighbor hopping integral and U characterizes 
the on-site interaction between particles. 

Eq.(pQ) actually can describe the behaviors both bosons and fermions, i.e., two spinless 
bosons or two electrons with opposite spins. 

For spinless bosons with the on-site interaction, the tight-binding Hamiltonian can gen- 



3 



erally be described by 

N-l N 
H = t Yl ( C n C "+l + C i+l C n) + E SnC n Cn 



n=l n=l 

N 



+ Uj2(clCnMc n ), (2) 



n=l 



where S n (c n ) is the boson creation(annihilation) operator of the nth site. The generic 
eigenstate for two spinless bosons is the superposition 



N N 

\(\] 

n\<ri2 ni<n2 



) = E Cna K"a> = E Cn 2 44 2 1°) > ( 3 ) 



where |0) is the vacuum and t/>" n2 is the amplitude of wave function. From Eqs. (J2J) and 
we can obtain the eigenvalue equation ([!]). 

For electrons with the on-site interaction, the tight-binding Hamiltonian can be described 

by 



N 



n=l a 

N N 

+ E E £ ^ C n,a C n,a + U E C i,T C ™,T C i,| C n,| , ( 4 ) 
n=l (j n=l 

where cJj CT (c nj(T ) is the electron creation(annihilation) operator for the electron at the nth 
site with spin o = ±|. For two electrons the Hilbert space can be conveniently divided 
into a singlet subspace with total spin 5 = and a triplet subspace with total spin S = 1, 
respectively. Since the triplet subspace permits no double occupation, it is not affected by 
the on-site interaction. In order to analyze the effect of the on-site interaction U, we will 
consider the case of the two electrons with opposite spins, i.e., the singlet subspace. In a 
chain with N sites the singlet subspace is spanned by iV(iV-|- l)/2 spatially symmetric basis 
functions 



J_(J J i J J 

t J * — - 



|ni,n 2 ) 

!.'-;,. tel. i + cj^cj^) |0) for ni < n 2 

( c L,T c La) 1°) for n i = n 2 

where |0) is the vacuum. An eigenstate for two electrons with the spatially symmetric wave 
functions is in general the superposition 

N 

I") = E ^i,«a K'^) > (6) 
ni<ri2 



where n2 is the amplitude of wave function. By making the transformation 

d,n 2 f O r n l < n 2, 



"•2 



^ V20Si,na fOT W l = ^2, 

the eigenvalue equation obtained from Eqs. (jU— [6]) can be written as equation 

In the following our numerical method is described for bosons. The extension for fermions 
is straightforward. 

B. von Neumann entropy 

For the two particles in the system we are studying, there are three local states at each 
site, |2) n , |l) n , |0) n , corresponding to the state with two, one or zero particles at the nth 
site, respectively. The local density matrix p n is defined Q,Q by 

Pn — Z2n |2) nn (2| + Z\ n | 1) nn (1| 

+ (l-zi„-z 2 n)|0>„„(0|. (7) 



For two spinless bosons, 



and 



Z2n = H c j n c n clc n I a) = Cn^C ( 8 ) 



Zln = (« | 4 C n I") - 2 ^2n 

JV N 

= E c, m <r m + E <>cv (9) 

m(m>n) m'(m'<n) 

Consequently, the corresponding von Neumann entropy related to nth site is 



E* n = -(1 - 2i n - Z 2n ) log 2 (l - Zi„ - ^2n) 

- Z ln log 2 2 ln - 2 2n log 2 Z 2n . (10) 

For nonuniform systems, the value of E° n depends on the site position n. At an eigenstate 
a, we define a site- averaged von Neumann entropy 



N 

E: 



i N 

v Ec en) 



N 

n=l 



For a delocalized state that all ib% „ = , 1 = for all m < no, this definition gives 
w -^log 2 y at large N, while for a localized state that ?/>"o n « = 1 at given n\ and n^, 
E£ = 0. In this paper, all the values of E% are scaled by j| log 2 y. From the two examples, 
we know that the scaled E® is close to 1 for eigenstates are extended and almost vanishes 
for eigenstates are localized. Henceforth, we omit "scaled" for simplicity. 

In order to analyze the influence of system parameters like the on-site interaction U, 
on the von Neumann entropy for all the eigenstates, we define a spectrum-averaged von 
Neumann entropy further gross measure 

a 

where M is the number of all the eigenstates. 



III. NUMERICAL RESULTS 



From now we consider only the repulsive interaction (U > 0). We directly diagonalize 
the eigenvalue Eq.flT]) with the periodic boundary condition at finite system sizes and obtain 
all eigenvalues E a and the corresponding eigenvectors \a). Without loss of generality, the 
hopping integral t is taken as units of energy. From the formulas (171— H2]L we then can 
obtain the site-averaged von Neumann entropy E" and the spectrum-averaged von Neumann 
entropy (E v ), respectively. 



A. TIP in a disordered potential chain 

For TIP in a disordered potential chain, the on-site potential e n in Eq.Q are random 
variables uniformly distributed with [— W, W\. Here W characterizes the degree of on-site 
disorder as in the Anderson model, j^] For this model, we calculate the spectrum-averaged 
von Neumann entropy (E v ) with a given set of parameters W and U. For every set of 
parameters W and U, the disorder average is taken over 100 samples. [271] More samples 
simply give similar results. 

For ID Anderson model in the absence of the interaction([/ = 0), it is well known that all 
the eigenstates are 



energy band center. 



ocalized and the one-particle localization length is £ ~ 25t 2 /W 2 at the 
FigiH gives the spectrum-averaged von Neumann entropy (E v ) versus 



6 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
W 

FIG. 1: The spectrum-averaged von Neumann entropy (E v ) versus the disorder W for 

U = at iV = 90. 



& u=o 




FIG. 2: The relation between the site-averaged von Neumann entropy E£ and the 
corresponding £ Q at W = 1 and U = for a typical sample. The E% versus log 10 is 

plotted in the inset. 

the disorder parameter W. It shows that (E v ) monotonically decreases as W increases, 
reflecting the trivial localization effect of the on-site disorder in the model. To find the 
correlation between (E v ) and the localization properties for TIP systems, we also study the 



A' 



inverse participation ratio (IPR)[l5|], defined by £ Q = ( N 2 Yl iV'nmJ 4 ) which gives 



the ratio of lattice sites occupied by particles to all lattice sites at an eigenstate a. The 
larger £ Q is, the more delocalized the eigenstate is. In FigfS] we plot the site-averaged von 
Neumann entropy E° versus £ a at U = for a typical sample with W = 1. On the whole, 
E® increases logarithmically with £ a as can been seen in the inset of FigJ2l so von Neumann 
entropy can well reflect the localization properties of two-particle eigenstates. 

For U > 0, the spectrum-averaged von Neunmann entropy (E v ) as a function of U at 
W = 1, 2, 3 are plotted in Fig 131 The results are similar for other W. For all W, we find that 



7 
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FIG. 3: The spectrum-averaged von Neumann entropy (E v ) as a function of interaction 
strengths U for {a)W = l,(b)W = 2 and (c)W — 3, respectively. Here N = 90. 

(E v ) first increases until its peak as U increases from zero, then decreases as U gets larger. 
There is a U*, at which (E v ) is equal to that at U = 0. When U is smaller (greater) than 
U*, (E v ) is larger (smaller) than the value of (E v ) at U = 0. To understand the effect of the 
interaction U on the von Neumann entropy, we calculate the site-averaged von Neumann 
entropy E° and illustrate the results in FigJH at different U for W = 1. It shows that 
only a small portion of the two-particle eigenstates have its value of E® changed by the 
interaction U. For small U, most of the newly created eigenstates are in the main band 
and the corresponding E° are larger than that in the noninteracting case, while for large U, 
most of them are above the top of the main band and the corresponding E" are smaller than 
that for particles without interaction. Therefore, the varying (E v ) with U is not monotonic 
as shown in FigJSJ These two competing effects reach an equilibrium at U — U*. 



For TIP in ID disordered potentials with U < t Shepelyansky [6j and others III, ll2] 
have found that the Anderson localization length increases with U . It agrees with our 
conclusions that for small U, (E v ) increases with U. For large interaction U, Evangelou et al. 
[lo| have pointed out that stronger localization occurs when compared to the noninteracting 
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FIG. 4: The site-averaged von Neumann entropy E£ and the corresponding eigenenergy 
E a at different U for a typical sample at W — 1. Here N = 90 and the spectrum-averaged 
von Neumann entropy = 0.8431, 0.8556 and 0.7916 for [7 = 0,2 and 6, respectively. 

case, because at U » t interaction can significantly modify the energy spectrum 6] and 



the created eigenstates have the comparatively small localization lengths. 



1 01 ] The result is 



consistent with ours that for U > U*, the values of (E v ) are smaller than that at U = 0. 



B. TIP in a quasiperiodic potential chain 



From Eq. ([TJ , the eigenvalue equation _ 
the Harper model can be described by [7 



or TIP in a quasiperiodic potential chain based on 



[A cos(27rcrni + Pi) + Acos(27rcm 2 + (3 2 ) + U5 niiri2 )ip nitn2 

= Eip nun2 , (13) 

here the parameter A characterizes the strength of the quasiperiodic potential, a and (3\^ 
are constants. As a typical case, we set a = (y/E — l)/2, (3i = (3 2 — 0. As is customary 
in the context of quasiperiodic system, the value of a may in fact be approximated by 
the ratio of successive Fibonacci numbers — F n = F n _2 + -Fn-i • in this way, choosing 
a = F n _i/F n m (\/5 — l)/2 and the system size = F n , we can obtain the periodic 
approximant for the quasiperiodic potential. In our calculation, N is chosen as Fibonacci 
numbers 34, 55 and 89, respectively. 

In the absence of the interaction U it is found 



28 



29| that that for A < 2 the spectrum 



becomes continues and all eigenstates are extended. For A > 2 the spectrum is pure point 




(a) * (b) * 

FIG. 5: (a)The spectrum-averaged von Neumann entropy (E v ) as a function of 
quasiperiodic potential strength A at U — 0. Here iV = 89. (b)The (E v ) versus A at U — 

for different iV. 

and all eigenstates are exponentially localized. For A = 2 the situation is critical with 
a singular-continuous multifractal spectrum and power law localized eigenstates. Metal- 
insulator transition can occur at A = 2. Obviously they are different from these for the ID 
Anderson model where all eigenstates are localized for W ^ 0. 

Figj5](a) shows the spectrum-averaged von Neumann entropy (E v ) with respect to the 
quasiperiodic potential strength A at U = 0. We observe that there is a sharp decrease 
of (E v ) near the critical value A c = 2. For A < A c which corresponds to extended states, 
all (E v ) are around 1, while for A > A c corresponding to localized states, all (E v ) are far 
less than 1. In Figfj^b) we plot (E v ) as a function of A at U = for N = 34, 55 and 89. 
The curves cross at A ~ A c . The crossing point separates the extended(A < A c ) and the 
localized(A > A c )regimes. 

For U > 0, the spectrum-averaged von Neumann entropy (E v ) as a function of U at 
different A are shown in FigJBJ Here we choose A = 1,2 and 3 as examples for extended, 
critical and localized regimes, respectively. At A = 1 and A = 2, (E v ) monotonically 
decreases when U increases from zero. At A = 3 (E v ) increases for small U and decreases 
for large U, respectively, which is similar to that shown in Figj3] for TIP in a disordered 
potential chain. For A in extended and localized regimes, (E v ) as a function of U are similar 
to that for A = 1 and A = 3, respectively. We also find that at the extended and critical 
regimes, the newly created eigenstates due to the interaction have small E®, which will 
reduce the value of (E v ). At the localized regime, the effect of interaction U on (E v ) is 
similar to that in the disordered potential system, so (E v ) increases for small U, decreases 
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FIG. 6: The spectrum-averaged von Neumann entropy (E v ) as a function of interaction 
strengths U for (a)A = l,(b)A = 2 and (c)A = 3, respectively. Here N = 89. 



16| found that the interaction 



for large U, respectively. 

For TIP in the ID Harper model, Shepelyansky et al. 
would induce localization effect for all A . At the same time, Evangelou and Katsanos 
js] found that, in the extended(A < 2) and the critical(A = 2) regimes, the velocity and 
the diffusion coefficient of TIP will decreases due to the localized pairing states. In the 
localization regime(A > 2), they found that propagation enhancement for small interaction 
and strong localization for large interaction, as in disorder systems. Apparently there are 
discrepancies in their results at the localization regime. Comparing to our results, in the 
extended and the critical regimes, we find (E v ) decreases with U, which agrees with both 
of the, concisions that the intern can induce Nation effect.Q |fa Q In the 
localized regime, we find (E v ) increases for small U and decreases for large U, respectively, 
which is consistent with the results of Evangelou and Katsanos [8] that the interaction has 
different effects on localization properties at small and large U. 
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FIG. 7: The spectrum-averaged von Neumann entropy (E v ) as a function of potential 
strength A at U = 0. The dashed line is linearly fitted for the corresponding data at 
< A < 2. Here na = 0.2, v = 0.7 and N = 89. 

C. TIP in a slowly varying potential chain 

Nest we study TIP moving in a ID system based on the slowly varying potential model 



24 



2a 



26I ] From Eq.flT]), the eigenvalue equation can be described by 



[A cos(7ran^ + Pi) + X cos(nari2 + (3 2 ) + U S ni;n2 ]ip ni>n2 

+^(^1+1,712 + V'ni — 1,712 + ^711,112+1 + V'ni.ns— l) 



(14) 



here A, a and v are positive numbers. For a irrational with -u > 2 or v = 1, this is 
equivalent to the models discussed in Sec. DT0 (A) and (B), respectively. For < v < 1, in 
the absence of the interaction U, it is well known [24j that there are two mobility edges at 
E c = ±(2.0 — A) provided A < 2.0. It is found that extended states are in the middle of the 
band ( \E\ < 2.0 — A ) and localized states are at the band edge ( 2.0 — A < \E\ < 2.0 + A 
). In other words, for A < 2.0, the extended and localized eigenstates coexist in contrast to 
the models studied in Sec. IIHI (A) and (B). For A > 2.0, all states are found to be localized. 
So in this model there are always localized eigenstates at A ^ 0. Obviously the spectrum 
properties are different from that for the Anderson model and the Harper model. 

Figd shows the spectrum-averaged von Neumann entropy (E v ) as a function of A at 
U — 0. For < A < 2 the data can be well fitted into a line. For A > 2 all data 
points lie far away from the proposed line, i.e., there is an abrupt decrease in (E v ) at the 
critical parameter A = 2. The critical parameter separates the localized regime(A > 2) 
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FIG. 8: The spectrum-averaged von Neumann entropy (E v ) as functions of interaction 
strengths U at different A. Here not = 0.2, v = 0.7 and iV = 89. 

and the regime where the localized and extended eigenstates coexist (A < 2). This result 
is consistent with the spectrum properties for the model, 24j so the spectrum-averaged von 
Neumann entropy (E v ) can describe the localization properties for two-particle systems. 

FigJH]shows the spectrum-averaged von Neumann entropy (E v ) versus the interaction U at 
A < 2. The results are similar to that for A > 2. We again find that for all A, as U increases, 
(E v ) increases for small U, decreases for large U, respectively, similar to that shown in Figj3] 
for TIP in a disordered potential chain. It indicates that propagation enhancement for small 
interaction and strong localization for large interaction. In other words, as long as there 
are localized eigenstate at single-particle models, the effect of interaction on the localization 
properties of TIP is similar to that in disorder systems. 



IV. CONCLUSIONS AND DISCUSSIONS 



With the help of von Neumann entropy, we have studied the effect of the on-site inter- 
action U on the localization properties of TIP in ID disordered, quasiperiodic and slowly 
varying potential systems, respectively. 

For TIP in a disordered potential chain and a slowly varying potential chain , as U 
increases from zero, we find that at first (E v ) increases for small U, then decreases as U 
gets large. It means that there are propagation enhancement for small interaction and 
stronge localization for large interaction. For TIP in ID Harper model with particles in the 
extended(A < 2) and critical(A = 2) regimes, (E v ) decreases as U increases, which indicates 
that the interaction would induce localization effect, while in the localized regime(A > 2), 



13 



the (E v ) as a function of U is similar to that for TIP in the disordered potential chain. From 
our studies, we find that (E v ) is a suitable quantity to describe the localization properties 
for two-particle systems. 

Summarizing all results from the three models, we can conclude that provided localized 
eigenstates at single-particle case exist, the delocaliation (localization) effect happens for 
small interactions (large interactions), while single particle states are extended or critical, 
the interaction always induces the localization effect. According to our results, we propose 
a consistent interpretation of the discrepancies between previous numerical results. 
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